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Abstract 

We consider the inhomogeneous generalization of the density matrix of 
a finite segment of length m of the antiferromagnetic Heisenberg chain. 
It is a function of the temperature T and the external magnetic field h, 
and further depends on m 'spectral parameters' ^j. For short segments of 
length 2 and 3 we decompose the known multiple integrals for the elements 
of the density matrix into finite sums over products of single integrals. This 
provides new numerically efficient expressions for the two-point functions 
of the infinite Heisenberg chain at short distances. It further leads us to 
conjecture an exponential formula for the density matrix involving only a 
double Cauchy-type integral in the exponent. We expect this formula to 
hold for arbitrary m and T but zero magnetic field. 
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1 Introduction 

The study of correlation functions of solvable quantum systems took a new direction in 
1992 when Jimbo et al. [16] shifted the focus of attention from the two-point functions 
to the density matrix of a sub-system as the object of principal interest. They managed 
to derive an m-fold integral expression for the density matrix of a chain segment of 
length m of the infinitely long, antiferromagnetic XXZ-spin-j chain in the off-critical 
regime. In the following years this work was generalized to the critical regime [17], 
to non-zero magnetic field [18] and, most recently [12, 14], to non-zero magnetic field 
and temperature. 

The above mentioned works have in common that they all deal with an inhomoge- 

neous generalization of the density matrix obtained by placing 'spectral parameters' or 
'inhomogeneities' onto m consecutive vertical lines of the corresponding vertex model. 
The homogeneous limit, when all the spectral parameters go to zero, is involved and is 
usually only taken at a late stage of the calculations after all formulae have been trans- 
formed into an appropriate form. Moreover, the inhomogeneous density matrix (at zero 
magnetic field and zero temperature) satisfies a first order difference equation in the in- 
homogeneities which is part of the novel reduced quantum Knizhnik-Zamolodchikov 
equation (rqKZ) [2,3]. 

In [5, 6] the multiple integrals for a specific zero temperature density matrix ele- 
ment of the isotropic Heisenberg chain were explicitly evaluated for m = 3 and m = 4. 
This line of research was further pursued in [7, 20]. In a second line of research repre- 
sentations of the zero temperature density matrix which do not involve multiple inte- 
grals were developed starting from functional equations [2,3,8-10,21,22]. This helped 
to work out further concrete examples of short-range correlation functions, such that 
nowadays closed analytic expressions for the ground state spin-spin correlators are 
known up to the seventh neighbour. In its most elaborate form of the rqKZ equation 
the functional equations were then used in [1] to derive an exponential formula for the 
density matrix of a segment of length m of the XXX chain that involves only a dou- 
ble Cauchy-type integral in the exponent. In the recent paper [4] such a formula was 
obtained for the XXZ and XYZ models as well. 

This article is concerned with an attempt to generalize both of the above approaches 
to finite temperature and partially also to finite magnetic field. We concentrate on the 
isotropic Heisenberg chain which on a periodic lattice of L sites has the Hamiltonian 

H = yt(oX+i-i)- (1) 

Here o", a = x^y^z, acts as a Pauli matrix on site j of the spin chain, and implicit 
summation over Greek indices is understood. For this model we have achieved a fac- 
torization of the density matrix for m = 2 and m = 3 into single integrals valid for 
arbitrary finite temperature and finite magnetic field. Based on this result and on re- 
sults obtained from the high temperature expansion of the multiple integrals [24] we 
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present a conjecture for a finite temperature generalization of the exponential formula 
in[l]. 

The paper is organized as follows. In section 2 we recall the multiple integral for- 
mula for the density matrix elements. In section 3 we explain the reduction procedure 
of those integrals in the cases m = 2 and m = 3. In section 4 we discuss our conjecture 
on the exponential form for zero magnetic field. In the appendix we show compact 
formulae for the density matrix elements for m = 3. 

2 The multiple integral formula for the density matrix 

The density matrix is a means to describe a sub-system as a part of a larger system in 
thermodynamic equilibrium in terms of the degrees of freedom of the sub-system. In 
our case the sub-system will consist of m consecutive sites of the spin chain. We first 
define a 'statistical operator' by 

PL = ^xp{-{H-hS')/T), (2) 

where = j'£!j=i^) is the conserved z-component of the total spin, T is the tem- 
perature and h is the homogeneous external magnetic field. In terms of this statistical 
operator the density matrix of a finite sub-chain of length m of the infinite chain is 
expressed as 

D(T\h)=\im^^^^^±^^. (3) 
L^'- tri...LpL 

By construction, the thermal average of every operator A acting non-trivially only on 
sites 1 to m can then be written as 

{A)T,h = tri...mAi,.,mD{T\h), (4) 

where Ai . is the restriction of A to the first m lattice sites. 

The tensor products ' ■■ ■ Sp'" composed of 2 x 2 matrices e^, a, P = 1 , 2, 
with a single non-zero entry at the intersection of row |3 and column a form a basis of 
End((C^)'^'") . The matrix elements of the density matrix with respect to this basis can 
be represented as 

SI viSm 

where the expression under the limit is the inhomogeneous density matrix element 
mentioned in the introduction. Due to the conservation of it is non-zero only if 
L7=i(a,-p,) = 0. 

For the non-zero inhomogeneous density matrix elements (5) of the XXZ chain a 
multiple integral formula was obtained in [12, 14]. When specialized to the isotropic 
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Figure 1 : The canonical contour C surrounds the real axis in counterclockwise manner 
inside the strip — ^ < ImX< j. 

limit (1) it reads 



n 



dOij 



11/ Yi(^j-^k-i) n i^j-^k) 



n 



dco^ 



detG(co;,^^) 



(6) 



The formula involves certain positive integers a+ and |3j derived from the sequences 
of indices (|3„)™^i specifying the matrix element. The indices take values 

1 , 2 corresponding to spin-up or spin-down. We shall denote the position of the j-th 
up-spin in (a„)™^j by aj" and that of the ^-th down-spin in {^n)n=i by • Then, 
by definition, = ^|lx+| j+i' where |a+| is the number of up-spins in (0Cn))f=i> and 



7 = 1, 



a 



+1 



Similarly (3^- =p^. |^+|, j= |a ' | + i, . . . ,m. 

The functions o(co), a(co) and G((0,^) are transcendental and are defined as solu- 
tions of integral equations. Through these functions and through the 'canonical con- 
tour' C (see figure 1) the temperature and the magnetic field enter the multiple integral 
formula. The integral equation for a(to) is non-linear. 



innru ^ , 2.7 /■ do) ln(l + a(co)) 

Ina(X) - -- + - - ^-p^^^-^ 



(7) 



and a((o) — 1/ a((o) by definition. Looking at the free fermion limit A — > of the XXZ 
chain [15] it is natural to interpret the combinations 1/(1 + a(to)) and 1/(1 + a(a))) 
as generalizations of the Fermi functions for holes and particles, respectively, to the 
interacting case. 
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The function G(to,^) is related to the magnetization density. It satisfies the Unear 
integral equation 



„ 1 _ r dO) 

^ '^^^ {X-Q{l-^-i) ~ Jc 7i(l + a(a))) 1 



)) l + (X-a))2 



(8) 



3 Reduction of multiple integrals 

The multiple integrals (6) can be reduced by means of the integral equation (8). We 
demonstrate the method with the example m — 2 and then merely present our results 
for m — 3. 



Emptiness formation probability for m = 2 

Consider (6) for m = 2 and ai = a2 = Pi = P2 = 1. Then 



Dil(^i,^2)fe-^i) = 



c 27i(l + a((Oi)) Jc 27i(l + a((02)) 



.det(G(co„^,)) (^i-ji-i)(^2-^2) ^ 



©1-0)2-1 



=:r(coi,(02) 



Here we may replace r{(Oi , (02) by [r((Oi , (O2) '"((02, C0i)]/2, since the domain of in- 
tegration of the double integral is symmetric. Now 



r(0L)i,0L)2) -r(co2,coi) = 

(0)1 -^i-i) ((02-^2) , ((02-^i-i)(wi-^2) ^((01,0)2) 



■ + 



CO1-CO2-I t0i-C02 + i l + (t0i-t02)^ 

and the polynomial P{(£)i , CO2) can be decomposed in such a way that 

/'(coi,co2) p(o)i)-;;(co2) 2 

--((01-0)2) 



, (10) 



with 



Then 



1 + (0)1-0)2)2 1 + (0)1-0)2)2 3 



i?(0)) = -0)3 - (^1 +^2 + i)o)2 + [i(^i +^2 + i) +2^l^2]0). 



(11) 



(12) 



Di}(^i,^2)(^2-^i: 



1 V- ■ M(0iG(0)i,^/>i) f d(i)2G{(02,^P2] 

- sign(P) 



c 7i(l + a(o)i)) Jc 7i(l + a(o)2)) 



p(O)i) 



1 + (0)1-0)2)2 3 



2 



(13) 



6 



The second term in the square brackets on the right is already of factorized form. The 
first term can be reduced to factorized form by means of the integral equation (8). 
Finally 



Z)}}(^i,^2)(^i-^2)= I sign(P) 



^ (3^Pi - ^ra + i)^i {^Pi ) + \fi2 {^Pi) 



where we have introduced a function 

¥(^1,^2) = 
and a family of 'moments' 



c7i(l + a((o)) (CO - ^2) (03- ^2 -i) 



^ r Jtoto-''~^G(to,E) 
^^^^^ = Jc .(l + a(co)) ' 



(14) 



(15) 



(16) 



Note that V(/(^i,^2) is symmetric. This can be shown by means of the integral 
equation (8). The physical meaning of \|/(^i, ^2) becomes evident in the limit of van- 
ishing temperature and magnetic field, where it can be expressed in terms of gamma 
functions. 



lim lim\|/(^i,^2) = 2idx\n 
r-^o/j— >o 



rG + f)r(i-f) 
r(i-^)r(i + f) 



(17) 



This is (up to a factor of —2) the two-spinon scattering phase [11,19]. It plays an im- 
portant role in the recent works [1,2] as it is the only transcendental function entering 
the general formula for the density matrix at zero temperature and zero magnetic field. 
Instead of y(^i , ^2) we shall rather use the closely related expression 



Y(^l,^2)=[l + (^l-^2fl¥(^l,^2)-l 



(18) 



in terms of which our final formulae look neater. We also define lim/,_»oY(^i)^2) =: 
Yo(^i,^2). 

Considering the moments in the same limit of zero temperature and magnetic 
field they turn into polynomials in ^ of order j —I, 



lim lim6y(E) = <bf\l) - (-idky~^ 



k=Q 



(19) 



for instance. 



,(0) 



r(^) = i, 4'^(^)=^+-, c^(^)=^^+i^. 



(20) 
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These polynomials satisfy the difference equation 

(^f (^-i)=2^^--i. (21) 

They allow us to define the 'nonnalized moments', 

(p/^)=(^x^)-^r(^)> (22) 

which vanish for T,h^O. We further introduce the symmetric combinations 

... ,._ det(cP;fe))|,,,^i,...,„ 

A„(^l,...,^„) — ^ , (23) 

lll<;<*;<nSA;; 

with the shorthand notation = ^fc — ^j- 

The A„ will turn out to be particularly convenient for expressing the density matrix 
elements for m = 2, 3. Using Ai and A2 in (14) we obtain 

} (^1 , ^2) = ^ + ^ (Ai (^1 ) + Ai (^2) ) + ^ A2 (^1 , ^2) - ^ Y(^i , ^2) . (24) 

The first moment cpi is exceptional among the cp^ in that it becomes trivial even for 
finite temperature if only the magnetic field vanishes. 



It follows that 



lim(pi(^) = 0. (25) 

/i— »o 



lim A^- (^) = , for all j e N. (26) 



Thus, for vanishing magnetic field, 

D}}(^i,^2) = ^-^Yo(^i,^2) (27) 
and, in the homogeneous limit ^1 , ^2 — 0, we have rederived the result 

(a5a|),,^„=4D||(7-|0)-l = i-limi / Z,.X,!:iTl (^^^ 

3 /!-»0 3 7t(l + Cl(tO)) tO(tO — 1) 

for the nearest-neighbour two-point function which alternatively can be obtained [13] 
by taking the derivative of the free energy with respect to 1/7. Here we can include 
the magnetic field into the calculation by simply taking the homogeneous limit in (24). 
We obtain 

{o\<f2)T,h = 4D\\{T\h)-2Ai{0) - 1 = ^A2(0,0) - ^y(0,0) (29) 

which seems to be a new result. Note that the magnetic field not only enters through 
A2 but also through y. The homogeneous limit of A„ exists for all n e N and is given 
by the formula 

lim . . . lim A„(^i , . . . , ^„) = det ^ . (30) 
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Complete density matrix for m = 2 

According to the rule tti + a2 = Pi + Pi which reflects the conservation of 5'^ the den- 
sity matrix for m = 2 has six non-vanishing elements. Using the Yang-Baxter algebra 
and certain identities of the type D\ (^) +D2(^) = 1 we find four independent relations 
between the six non-vanishing elements, 

d}2(^i,^2)=d}(^i)-D}}(^i,^2), Dl\i^,,^2)=D\{^2)-D\\{^r,^2), 

Z)li(^l,^2)-/>il(^l,^2) = ''^^^^^~''^^^^^ (31) 

1^12 

Thus, we have to calculate one more independent linear combination, say (^1 j ^2) + 
(^1 , ^2) 5 in order to determine the complete density matrix for m = 2. The calcula- 
tion can be done along the lines described above. 

In order to obtain a convenient description of all density matrix elements we shall 
resort to a notation that we borrowed from [4]'''. We arrange them into a column vector 
hm e (C^)®^"* with coordinates labeled by +, — instead of 1,2 according to the rule, 

^ej,...,s„,e„,...,s. (X,, . . . , = D^^;:::^^{^^,. n(-ey) , (32) 

where A,y = — i^j for j = 1 , . . . , m. 

Then, setting = (/iJ+"",/iJ"+",/iJ""+,/i2 ++",/i2 ^"^,^2 "^^)' 

1 1 
v= ^vo-— Y(^i,^2)vi 

+ i(Ai(^i)+Aife))v2-^(Ai(^i)-Aife))v3 + ^A2(^i,^2)v4 (33) 

with 

vo = (l,-l, 0,0,-1, 1)S 
Vl = (1, 1,-2, -2,1, 1)S 
V2 = (l,0,0,0,0,-iy, 

V3 = (0,l,i^r2''-i^r2:-l:0)': 

V4 = (l,l,l,l,l,l)^ (34) 

From this we can read off the transverse neighbour correlation functions in a magnetic 
field to be {&{G'2)T,h = -5A2(0,0) - ^7(0,0). 

'^This definitioii was first introduced in [2] and later modified in [4]. 
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Emptiness formation probability for m = 3 

As in the case m = 2we can reduce the triple integral (6) representing the emptiness 
formation probability for m = 3 to sums over products of single integrals. Again we 
have to decompose the rational functions in the integrand appropriately and use the 
integral equation (8). It turns out that the final result can be represented in terms of the 
functions y and Aj, j — 1,2,3, 

i'lllSl.fc,^3) = 

24 40^12^13 24^13^23 60 

. 1-^13|23 ,e e X 3 + 2^^2 +^^13^23 \\ 1^ \ 

24^13^23 120^13^23 



+ cyclic permutations. (35) 
In the limit of vanishing magnetic field (26) applies and our result reduces to 



D}i}(^i,^2,^3) = :^ + ^r-^^Yo(^i,^2)+cyclic permutations. (36) 

24 24^13^23 

Note that the only effect of taking the limit T ^ is that the function yo (^i , ^2) changes 
into its zero temperature form (17), (18). 



Complete density matrix for m = 3 

We have factorized the full density matrix for m = 3. It contains 20 non- vanishing 
matrix elements, each of similar form to that in (35). Except for the emptiness for- 
mation probability Z)||| we also reduced the integrals for the symmetric combination 
^121+^211 +^121 +^211 +^112 +^112- Using relations similar to (31) and the high 
temperature expansion data for the inhomogeneous density matrix elements up to the 
order this was then enough to conjecture the complete density matrix for m = 3. 
We show it in a compact notation in the appendix. For our purposes here it is suffi- 
cient to know that all density matrix elements are linear combinations of the functions 
Ai, A2, A3, Y and yAi with coefficients rational in the differences '^jk- For the case of 
vanishing magnetic field an 'exponential formula' for all density matrix elements is 
suggested in the next section. 



The next-to-nearest neighbour two-point functions 

It is a rather straightforward exercise to work out the next-to-nearest neighbour two- 
point functions from our general result for the m = 3 density matrix. We have to take 
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the appropriate linear combinations of density matrix elements and have to carry out 
the homogeneous limit , ^2, ^3 — 0. We obtain, for instance, 

(olo|)r,/z = ^A2(0,0)-^Y(0,0) 

- ^(A2)x.(0,0) + ^(A2)xy(0,0) - ^Yxc(0,0) + ^Yxy(0,0) , (37a) 

Kc^)r,/,= -^A2(0,0)-iY(0,0) 

+ ^(A2)«(0,0) - ^(A2)xy(0,0) - ^Y«(0,0) + ^Yxy(0,0) . (37b) 

Here we denoted derivatives with respect to the first and second argument, respectively, 
by subscripts x and y. Equations (37) generalize an important result of Takahashi [23] 
to include the temperature and the magnetic field. 



4 The exponential formula 

We observed in the previous section that for m = 2, 3 the density matrix for zero mag- 
netic field is determined by a single transcendental function Yo (recall that lim/,_»o Ay = 
0). The situation is the same as for zero temperature. In fact, even the coefficients 
agree. Hence, it is tempting to substitute the function Yo for its zero temperature ana- 
logue into the general exponential formula recently obtained in [1]. This formula then 
gives the correct result for the 6 non-trivial density matrix elements for m = 2 and also 
for the emptiness formation probability and for the symmetric combination of density 
matrix elements mentioned in the previous section for m = 3. It further coincides up to 
order with the high-temperature expansion data for all 20 non-vanishing inhomo- 
geneous density matrix elements for m = 3. For m = 4 we compared the conjectured 
form in the homogeneous case with the high-temperature expansion obtained from the 
homogeneous version [12, 14] of the multiple integral formula (6). For the emptiness 
formation probability we found full agreement up to the order of r~^^ 

Conjecture 1. The density matrix of a finite sub-chain of length m of the infinite XXX 
Heisenberg chain at finite T (for h = 0) is determined by the vector 

hmQ^li • • • ■I'^m) — "' ' '" — \^\^ ^jj 1 (38) 



"-(^^'•••'^"'^ = ^^// 2^2^ [1-(a'i-A'2)2]2 ^^^^ 



X 



through (32). By the integral over n\, H2 it is meant to take the residues at the poles 
A,i , . . . , A,^ of the integrand. 
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For the notation we are referring to [1]": The vector s = (g) <8) (^) — (j) (8) (q) is 
the spin singlet in ® C^. The vector spaces in (C^)®^'" are numbered in the order 
l,2,...,n,n,n — 1,...,1. This defines s^. ~ is the projector onto the one-dimensional 
subspace of spanned by s. 

In order to define the transfer matrices in the integrand in (39) we first of all intro- 
duce an L-matrix L(k) e U (3(2) ® EndC^, 

L(^) = |^(2X+l+L«®a«), (40) 

where the E*^ G SI2 are a basis satisfying [E", E^] = 2i£'^P^E^, where d is determined by 
the Casimir element through d-^ = (E")^ + 1 and where p^k^d) satisfies the functional 
relation 

2-2'k-d 

p{X,d)p{X-l,d)= (41) 
(for more details see [1]). Then, for integer z, the 'transfer matrices' 

trj.r(A,;A,i,...,A,n) = 

tr.Lj (X - - 1 ) . . . (X - ?i„ - 1 )L„ (X - X„) . . . Li (X - ) (42) 

entering (39) are defined by substituting the irreducible representation of U{sl2) of 
dimension z into the definition (40) of the L-matrices. For non-integer z this can be 
analytically continued into the complex plane. 



5 Discussion 

Starting from the multiple-integral formula (6) we have investigated the density matrix 
of a finite segment of length m of the infinite isotropic Heisenberg chain at finite tem- 
perature and finite magnetic field. We found that the multiple integrals can be reduced 
to sums over products of single integrals, in much the same way as for T,h = 0. On 
the one hand this gives new efficient formulae for the calculation of finite-temperature 
short-range correlations of the XXX chain in the thermodynamic limit. On the other 
hand this shows that the density matrix and the correlation functions of the Heisen- 
berg chain at finite temperature and finite magnetic field may be explored in a similar 
manner and to much the same extent as in the ground state case without magnetic field. 
When the magnetic field is switched off, but the temperature is kept finite the 'algebraic 



llln fact, the only difference between our formula (38), (39) and the result of [1] is in the function Yq. 
In [1] a function (O was used which is related to Yo by 



T^0 2{\-{Xl-'k2f)' 
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structure' of the density matrix, as it shows up in the rational functions in the differ- 
ences of the spectral parameters, seems to be still the same as for zero temperature, and 
only the 'physical part' of the expressions, encoded in the transcendental function yo, 
changes. This is certainly true for m < 2 and for some of the density matrix elements 
for m = 3. From our high-temperature analysis it seems most likely true also for all 
density matrix elements for m = 3,4, whence our conjecture 1 in the previous section. 
A second conjecture we are tempted to formulate inspecting our results for small m is 
the following: In the case of non-vanishing magnetic field the density matrix elements 
for a segment of length m seem to depend only on y and on Ai , . . . , A„j. 
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Appendix: Density matrix elements for m = 3 

Here we show our result for other elements of the density matrix in the m = 3 case. 
We consider the 3 by 3 block 

/7)112 r)112 r)112\ 
/ ^112 ^121 ^211 \ 

D3x3= \D\l\ D\l\ D\l\ . (A.1) 
\r)2ll r)2ii n2ll / 
\'^112 -^121 ''^211/ 

Another 3 by 3 block can be obtained from (A. 1) by the substitution 1 <-> 2 and h — > —h. 
The answer looks as follows 

l<i<j<3 \<i<i<y,k^i,i 

+ j^BiAl{^i)+ I fiyA2(^,-,y+5i,2,3A3(^i,^2,^3), (A.2) 

i=l l<!<i<3 

where jij = y{^,i,^j) and Ao,A,- j,A,- y|^,5,-,5,-,;,5i,2,3 are 3 by 3 matrices. In order to 
define them let us introduce 3 by 3 matrices with the following elements 

-1 i{k^l = i 

{eo)k,i = {Zo)k,i = 1, {Xi)k,i = { 1 iik = l^i, 

otherwise 

X±. = e/±4. fori<;, (e/)^,/ = 5,-^5,-/, l</,;,A:,/<3. (A.3) 
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We also define Xj*] = xf^^ and X;+3„ = Z,, X^^^ j_^2m ~ ^tj for «,m e Z. Then 



^— ^1,2 ~ ^1,3 +^2,3' 



1 1 i _ 1 i _ 11/1 +\ 



4i(^-V.-(-i)'^--w-(-')'^-V.) 

Sl.2,3=-^Zo. (A.5) 

It is implied in the above formulae that the triple (/, 7, k) is always the cyclic permuta- 
tion of (1,2,3) and also that A,-,j = Aj ,-,A,j|^ = Aj^i]^j^,Bij = Bj i. Let us also mention 
that the matrices Aq and A, y coincide with those which can be obtained from the for- 
mula (39). 
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